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ABSTRACT 


Intensity  interferometry,  as  developed  by  Hanbury  Brown  and  Twiss 
for  stellar  observation,  has  shown  relative  .insensitivit)  to  atmospheric 
scintillation.  However,  with  clasiical  sources,  the  limitations  placed 
on  this  technique  by  quantum  noise  and  defector  efficiency  are  severe. 
This  situation  is  vastly  improved  when  lase  •  illumination  is  employed. 
Generalizing  a  form  of  the  mutual  coherence  function,  we  derive  the  far- 
zone  behavior  of  the  mutual  intensity  function  for  an  intermediate  time 
average.  This  result  is  used  to  reconstruct  the  i rradiar.ee  distribution 
of  a  spatially  rough  source.  The  far- field  intensity  distribution  is 
recorded  spatially  for  one  time- resolution  unit  of  the  detector.  The 
resulting  spatial  signal  is  auto  correlated  and  related  to  the  intensity- 
distribution  over  the  source.  Thus,  without  averaging  in  the  time 
domain,  a  spatial  Fourier- transform  relation  is  derived  between  the  far- 
field  intensity  correlation  and  the  source  irradiance,  similar  to  the 
results  of  Hanbury  Brown  and  TWiss. 
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I.  INTRODUCTION 


In  ar.  earlier  BRL  Report  under  this  title,  we  developed  the 
mathematical  relationship  between  t?  j  intensity  distribution  over  a 
rough  surface  and  the  intensity  correlations  in  the  far  field.  Although 
intrinsic  to  the  underlying  mathematics,  the  assumption  of  gaussian 
statistics  was  not  explicitly  used.  In  an  effort  to  extend  the  approach 
to  sources  of  arbitrary  degree  of  spatial  coherence,  methods  of  gaussian 
statistics  have  been  applied  to  this  study.  This  leads  to  a  simpler 
mathematical  development. 


The  final  results  of  this  report  are  nearly  identical  to  the 
* 

earlier  work.  However,  here  the  detailed  development  of  the  work  of 
Msrchsnd  and  Wolf  is  omitted.  We  have  retained  the  original  introduction 
so  thst  this  report  may  be  self  complete. 

■>b=  history  of  intensity  interferometry  is  rooted  in  the  work  of 
.'{anbury  !>own  and  Twiss ,  Their  earliest  investigations*  dealt  with 
the  problem  of  resolving  stellar  radio  sources  by  a  technique  involving 
The  correlation  of  the  squared  outputs  of  two  receivers,  The  advantages 
are  the  reduction  of  certain  kinds  of  experimental  constraints  as  well 
as  the  comparative  insensitivity  of  the  method  to  atmospheric  scintil- 
lation.~  A  preliminary  conclusion  reached  at  that  time  was  that  this 
technique  of  intensity  correlation  would  not  be  applicable  at  optical 
frequencies  because  of  limitations  imposed  by  photon  noise.  However, 

In  later  work,'  Banbury  Brown  and  Twiss  showed  that  meaningful  inten¬ 
sity  correlations  could  be  made  at  optical  frequencies  even  with  highly 
degenerate  sources.  The  limitations  placed  on  this  approach  by  quantum 
noise  and  detector  efficiency  have  been  severe,  calling  for  highly 
refined  experimental  technique. 


»  * 

'  ERL  Report  No.  1616 
** 

References  may  be  found  on  page  ??. 
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For  laser  illumination,  the  situation  is  very  different.  The 

signal-to-noise  ratio  can  be  typically  increased  by  six  orders  of 
4 

magnitude.  However,  the  statistics  of  the  source  must  be  considered 
in  the  measurement.  The  key  to  relating  intensity  correlations  to  some 
property  involving  field  correlation*  lies  in  the  assumption  of  gaussian 
statistics,^  for  which  all  higher  moments  are  determined  from  the  first 
and  second.  Single-mode  lasers,  though,  are  distinctly  nongaussian  in 
their  statistics,  and,  therefore,  cannot  be  described  by  theory  framed 
for  thermal  sources.  But  with  the  addition  of  only  a  few  axial  modes, 
the  field  amplitude  becomes  nearly  gaussian  distributed.^ 

The  principal  formula  used  by  Hanbury  Brown  and  Twiss^  to  infer 
the  diameter  of  a  distant  source  shows  that  the  time-averaged  corre¬ 
lation  of  intensities  at  two  points  is  equal  to  the  product  of  a 
function  involving  the  temporal  characteristics  of  the  source  with  the 
square  of  the  spatial-Fourier  transform  of  the  source  intensity  distri¬ 
bution.  If  the  intent  of  an  intensity  correlation  experiment  is  to 
gain  information  concerning  the  source  intensity  distribution,  then  the 
source  temporal  statistics  may  be  of  little  interest  in  themselves. 
Consideration  of  the  temporal  statistics  is  necessary  if  the  intensity- 
product  output  of  the  detectors  is  averaged  in  the  time  domain  (as  it 
nearly  always  is]  to  overcome  the  limitations  imposed  by  photon  and 
detector  noise  and  possibly  reduce  scintillation  effects  produced  by 
transmission  through  the  atmosphere. 

Intensity  interferometry  can  be  understood  as  a  two-point  corre¬ 
lation  of  intensities  following  the  squaring  of  the  electric  field  at 
the  detector.  If  the  source  is  quasi -monochromatic,  each  differential 
element  on  the  object  emits  a  number  of  temporal  modes  that  interfere 

g 

with  ether  at  the  detector.  If  the  detector  has  sufficient  speed 

to  detect  these  beat  frequencies,  the  amplitude  and  phase  of  the  in¬ 
coming  intensity  beats  are  utilizer.',  ft  is  commonly  argued  that  in 
these  circumstances,  the  random  flui.;uaticns  of  the  temporal  statistics 
from  different  points  of  the  source  cause  the  beat  frequencies  from 


S 


each  source  point  to  add  incoherently  at  the  detector.  "Hie  time- 

averaged  intensity  correlation  is  then  proportional  to  the  squared 

spatial-Fourier  transform  of  the  source  intensity.  This  approach  gives 

9 

essentially  a  squared  version  of  the  van  Cittert-Zemike  theorem. 

We  suggest  that  the  requirement  of  surface  roughness  at  the  source 
(to  assure  spatial  incoherence)  is  sufficient  to  guarantee  the  in¬ 
coherent  addition  of  beat  frequencies  at  the  intensity  detectors.  Thus, 
if  temporal  noise  (photon  noise,  time- dependent  detector  noise)  is 
largely  absent  in  a  local  spatial  sense,  as  might  be  the  case  with  a 
multi-axial-mode  laser  with  photographic  detection,  then  the  intensity 
information  might  be  gathered  during  one  resolution  time  of  the  detec¬ 
tor  over  a  plane  section  normal  to  the  direction  of  light  propagation. 
Any  noise  arising  in  the  process  would  be  spatial,  and  might  be  averaged 
out  by  taking  a  sufficiently  large  area  of  spatial  correlation.  The 
reduction  of  atmospheric  spatial  noise  would  be  similar  to  a  process 
known  as  aperture  averaging. 10  Film- grain  noise  would  be  extremely 
well  averaged  by  the  relatively  large  area  of  averaging. 

The  relative  insensitivity  of  intensity  interferometry  to  turbu- 
-  lence  is  due  to  the  assumed  dispersion less  nature  of  the  propagation 
medium.  Because  each  temporal  frequency  sees  the  same  refractive  index, 
the  differential  (beat)  frequencies  remain  unchanged.  However,  the 
spatial-Fourier-transform  relation  between  the  source  and  the  far  field 
is  scaled  as  the  average  frequency,  not  the  beat  frequency,  and  thus 
the  resolution  afforded  by  optical  frequencies  is  maintained. 

The  idea  of  examining  spatial  beat  frequencies  of  second-order 
correlation  is,  of  course,  not  new.  M_ny  classical  fie la- correlation 
interferferomeLers,  as  well  as  holographic  experiments  are  built  on  this 
principle,  involving  a  spatial  or  time  lag  between  interfering  beams  of 
the  same  source.  More  difficult  is  the  spatial  recording  of  beats  from 
two  independent  sources,  as  demonstrated  by  Magyar  and  Mandel. 
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We  have  examined  the  problem  of  relating  far-field  spatial  inten¬ 
sity  (fourth-order-field)  correlations  to  the  intensity  distribution  of 
the  source,  without,  here  considering  the  limitations  due  to  noise.  Our 
approach  has  been  completely  classical,  drawing  on  a  straight-forward 
generalization  of  a  method  given  recently  by  Marchand  and  Wo1 f .  ^  Our 
notation  is  similar  through  the  development  of  their  Eq.  (35). 


II.  THE  INTERMEDIATE-AVERAGE  MUTUAL  COHERENCE  FUNCTION 


For  a  stationary  scalar  wave  field,  the  mutual  coherence  function 
for  the  correlation  of  two  space-time  points  is  often  written 


r  Cx.2 2 » T  »T^ 


T 

VT(x1,t+t)  V*(x2,t)  dt, 
-T 


1 

2T 


(1) 


where  x  =  x  i  +  y  i,  and  the  limits  for  the  time  integration  are 
allowed  to  approach  infinity.  For  this  case,  however,  we  wish  to  keep 
the  parameter  T  finite,  and  by  the  subscripts  indicate  that  we  assume 
a  knowledge  of  VT(x^,t+-r)  and  VT*(x_2,t)  only  over  the  finite  sample 
length  2T.  We  wish  to  call  T  (Xj,x_2,f,T)  the  intermediate-average  mutual 
coherence  function  and  stress  that,  for  arbitrary  T  or  shift  of  origin, 
it  may  bear  little  resemblance  to  the  mutual  coherence  function  defined 
by  the  ensemble  average. 


Following  Ref.  i2,  we  represent  V(x,t)  as  the  temporal  Fourier 
transform  of  the  complex  analytic  signal  (where  a  factor  of  (2W)  *  has 
been  suppressed) , 


and 


VT(x„t) 


j>0O 

vT(x,0))  exp(-iu.'t)d{.>  for  0  <  (t{  <T 

0 


=  0  otherwise 


T 

vT(x,w)  =  I  VT(x,t)  exp(iut)  dt. 
-T 


(2a) 

(2b) 


(2c) 


10 


Substituting  Eq.  (2 a)  into  Eq.  (]V  interchanging  the  order  of  integra¬ 
tion  and  time-averaging,  and  performing  the  time  integration,  wo  get 


vj 

► 

T  j  >  1)  —  hy,  2 j  x2  > 


expf-itt^T)  sine  [(a)^ -(1)2 )T]  da^du^* 


where 


_  sm  x 

sine  x  -  - 

x 


and  the  function 


M-r-2  E  vT^-r“i-)  vt*  Qb'V 


is  termed  the  cross-spectral  density.  The  subscript  T  here  and  later 
implies  a  function  based  on  the  electric-field  statistics  only  for  the 
particular  sample  2T  in  length  about  the  origin.  The  sine  function  of 
Eq.  (3)  assumes  the  role  of  a  low-pass  filter.  If  T  is  very  small,  the 
two  frequency  variables  of  Eq.  (3)  are  essentially  independent  and  all 
cross  terms  are  represented  in  the  product  of  Eq.  (5).  These  cross 
terms  form  a  high-frequency  spectral  content.  However,  as  T  tends  to 
infinity,  the  sine  function  assumes  the  role  of  a  6  function,  constrain¬ 
ing  correlation  to  occur  only  between  identical  frequencies  in  the 
transform  product  and  forcing  the  integral  to  a  one-dimensional  form. 

In  the  limit  of  large  T,  the  filtered  spectrum  of  Eq.  (5)  becomes  the 
mean-square  value  (dc)  of  each  temporal -frequency  component  in  the 
signal. 

Omitting  here  the  details,  the  far  zone  form  of  the  cross-spectral 
density  function  nas  been  derived  by  the  method  of  Marchand  and  Wolf 
wi’.h  the  exception  that  the  parameter  T  is  left  finite.  The  resulting 
.'-’-.•im  of  the  intermediate-average  mutual  coherence  function  is 


r(£.p*2T,T-) 

co 

=  cos61  cos62  JJ  expC-iosjx)  sinc^u^-w^T] 


X  k,k  d^da^ 


exp[i(k1r1-k2r2)] 


rlr2 


/  X1  yl  X2  y2  \ 

X  VT  \^kl  rj*  kl  rx;  °;  k2  r2J  k2  r2;  0;  V^j  * 


(6) 


where  the  spatial- frequency  correlation  function  in  the  plane  z  =  0  is 
defined 

.  /  xi  >'i  \  *  *  /  x2  y2  \ 

vT^kl  7?  kl  ^  0;  “;)  VT  (k2  T2‘  k2  0;  "2j 

/  X1  yi  X2  y2  \ 

=  v„  k  -i,  k  r1;  0;  k_  k,  0;  o  ,  (7) 

‘  \  1  rl  1  rl  2  r2  2  r2  1  -) 


and  v^(f,g;z;w)  is  related  by  a  two-dimensional  spatial-Fourier  transform 
to  the  time- frequency  representation  of  the  electric  field  by 


vT(x,y,z;oj)  = 


vT(f,g;z;u)  exp [i (fx+gy) ]  dfdg. 


(8) 


Finally  the  circular  functions  give  the  direction  cosines  for  the  two 
correlated  fields,  while  the  arguments  are  the  angles  that  the  two 
directions  make  with  the  z  axis;  specifically 


zl  Z2 

—  =  cose,  and  —  =  cose,, 

Tj  1  r2  2 


(9) 
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If  T  is  allowed  to  approach  infinity,  Eq.  (6)  reduces  to  the  form  given 
in  Ref.  12  using  Eqs.  (5)  and  (35).  We  are  now  ready  to  form  the  self¬ 
intensity  function  in  the  far  field. 


III.  THE  SELF- INTENSITY  FUNCTION  IN  THE  FAR  FIELD 

We  now  examine  the  form  of  the  self-intensity  in  the  far  field  by 

letting  points  x^  =  =  —  ^d  then  lettin2  the  time  delay,  t,  be  zero. 

Under  these  conditions  the  mutual  coherence  function  reduces  to  the 
]  3 

self-intensity,  and  using  Eq.  (6),  we  have 


I(x,T)  =  cos^0 


[{ 

d^2^r 

i  c 


exp[i(k1-k2)r] 


sine  [Cw1-w2)T] 


x  v. 


(kl  f-  kl  r’  0;  “l)  V  (k2  f-  k2  r’  0;  “2)' 


(10) 


where,  as  indicated  earlier,  the  sine  function  acts  to  suppress  temporal 
frequencies  in  the  cross  spectrum  higher  than  *vl/(2T)  Hz.  We  now 
utilize  the  linear  transformation  of  the  time-frequency  variables  (for 
which  the  jacobian  is  unity)  defined  by 

u,'i+u>2 

Uj  -  is2  :  ?  and  — 2 —  =  o  •  (H) 


Writing  the  w  variables  in  terms  of  these  center-of-mass  coordinates, 
we  get 


u 

1 


2a+p 

2 


and  io? 


2  g-p 
2  * 


(12) 


which,  when  substituted  into  Eq.  (10),  gives 
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f 


x  cos  6 


sinc(pT)  exp(ipr/c)d 


vT  f .  o;  -p/2) 

»  v(^f-  Oio-p/z)  , 


where  the  dependence  of  the  amplitude  on  the  difference -frequency  co¬ 
ordinate,  p,  has  been  dropped,  since  for  quasi-monochromatic  radiation 

o2  »  |p2/4|. 

Now,  using  the  defining  transform  relation  of  Eq.  (8),  we  write 
the  spatial-correlation  function  at  the  source  (z  =  0)  where 


2i£/A  o:  0*0/2) 

\  c  r  c  r  / 

'  *  (a-o/2  x  o-p/2  y 
x  VT  r-  c  r  0! 


— ^—7  f  f  vTU.,n,;0;  o+p/2) 

(2ity  JJJJ  1  A 


vT’  (57,n2i°J  o-p/2) 
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exp 


5)] 


x  exp  1 


^)]  «lW„; 


(14) 


Now  the  intermediate-average  spatial-correlation  function, 

vT(51>n1;0;  i+ p/2)  vT*(c2,n7;0;  o-p/2)  , 

when  considered  wit!  che  filtering  action  of  the  sine  function  of  Eq.  (13), 
will  have  an  effective  contribution  only  for  the  low-frequency  components 
formed  by  the  difference- frequency  terms  -u  1/(2T)  Hz  or  less. 

In  addition,  we  assume  the  mode  population  to  be  a  slowly  varying 
function  of  o,  since  o  »  p/2,  and  thus  we  write 

[VT(C1’V0;  °+p/2)  VT*(W0;  °"p/2^low  freq. 


=  o+p/2)  A(s2,n2;  o-p/2) 


exp[i$  (Cj-?2>  r'1~ri2;p) 


(15a) 


where 


=  AUj.o)  AU2,o)  H(p)  exp[i$  (j^-j^p)] 


H(p)  = 


1  for  A(C;o±p/2)  ?  0 
0  otherwise. 


(15b) 
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We  note  that,  in  Eq.  (lSb) ,  a  spatial-vector  symbolism  has  been  adopted 
to  shorten  the  notation,  where  +  pj.  Equation  (15)  indicates 

the  loss  of  the  optical-frequency  phase,  whereas  the  phase  of  the  in¬ 
tensity  envelope  formed  by  temporal  beat  modes  is  preserved.  The  degree 
to  which  the  phase  of  this  envelope  is  detected  depends  on  the  bandwidth 
of  the  source  and  the  detector  resolution,  2T. 

Essentially,  the  approximation  of  Eq.  (15b)  was  made  by  lianbury 
Brown  and  Twiss,  except  for  the  defining  of  the  H  function.  Its  intro¬ 
duction  is  brought  about  by  the  description  of  narrow-band  sources  by 
terms  in  A(^,o).  For  a  thermal  source  of  relatively  large  bandwidth, 
many  of  the  higher  beat-frequency  (p)  components  will  not  be  resolved 
by  the  system  response  (s,2T)  described  by  the  sine  function.  In  addi¬ 
tion,  the  beat  spectrum  will  be  a  continuous  function.  For  a  laser 
source  exhibiting  a  series  of  axial  modes,  however,  the  complete 
difference- frequency  domain  might  lie  entirely  within  the  system  re- 
snonse  but  be  piece-wise  continuous  in  its  extent. 

Relative  to  tha  representation  of  the  intermediate  average  by  'he 
form  of  Eq.  (15),  we  wish  to  reiterate  a  statement  made  following  Eq.(l) 
that  the  intermediate-averaging  process  may  bear  little  resemblance  to 
the  infinite  time  average,  even  so  far  as  the  detail  of  the  ampl*‘  a 
terms,  A(,£,a).  fh*s  situation  would  be  serious  if  our  intent  were  to 
infer,  fo7  example,  the  time -frequency  statistics  of  the  source.  But 
in  the  present  concept,  we  desire  only  to  infer  the  spatial  properties 
of  the  source.  If  we  consider  s  multi-axial-mode  laser  beam  scattered 
from  a  spatially  rough  surface,  the  lack  of  correspondence  between  the 
two  averages  is  unimportant,  for  all  such  mode  history  is  integrated 
out;  all  areas  of  the  scatterer  see  the  same  mode  characteristics.  Any 
mode  fluctuation  would  be  seen  simply  as  a  variation  in  total  received 
power  from  one  sample  to  the  next.  Here,  we  simply  require  for  one 
detector- resolution  time  over  a  spatial  domain  that  the  process  of 
Eq.  (15)  exhibit  a  minimum  of  two  temporal  modes  [to  maintain  the  phase 
term  ^(^j-j^P)  with  sufficient  mode  population  (reflected  in  the 
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amplitude  terms  A(£,o) }]  such  that  quantum  noise  in  both  the  carrier 
wave  and  the  detector  can  be  ignored. 

Using  the  results  of  Eq.  (15)  in  Eq.  (13)  and  taking  0  «  1,  we 

write 

OO  CD 

# 

c"da  dp  exp(ipr/c)  sinc(pT)  H(p) 

0 

CO 

*  ||  ACI^0)  A (12,0)  expriiKij-j^Jp)] 

-CD 

exp  dl'i)] 

exp  [i(2;fr-)  (h-$  diid^  •  (I6) 

We  now  have  the  self-intensity  in  the  far  field  as  a  double  integral 
over  sum-  and  difference -frequency  components  as  well  as  two,  two- 
dimensional  spatial  Fourier  transforms  over  the  source. 

14 

Following  Goodman,  we  argue  that  the  received  field  at  any  point 
in  the  far  zone  consists  of  a  sum  of  random- amplitude,  random-phc°' , 
complex  phasors  contributed  by  the  elementary  scatterers.  If  the  size 
of  the  scattering  area  is  large  enough  to  include  many  point  scatterers 
(or  there  are  enough  elementary  coherence  areas  composing  the  source), 
the  central  limit  theorem  may  be  used  to  conclude  that  the  electric 
field  in  the  detection  plane  is  <■  gaussian  random  piocess  in  a  spatial 
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sense. 


Using  the  form  of  Eq.  (lb)  and  its  property  of  spatial  gaussian 
statistics,  we  are  ready  to  form  the  fourth-order  correlation  function 
in  the  far  zone. 

IV.  FOURTH -ORDER  FIELD  CORRELATION  IN  iUE  FAR  ZONE 

Basic  to  the  theory'  of  intensity  correlation  is  the  ability  to 
represent  higher- order  field  correlations  (fourth,  sixth,  etc.)  in 
terms  of  tie  first  and  second.  This  is  due  to  a  well-known  property  cf 
gaussian  processes.  However,  in  the  work  of  Danbury  Brown  and  Twiss, 

this  property  was  never  explicitly  used,  although  it  was  discussed  later 
by  Kolr  as  a  plausibility  argument  for  their  work. 

To  form  the  fourth-order  correlation  function,  wo  can  proceed  1))' 
writing  the  two-point  product  of  intensities  in  the  far  field  using 
Eq.  (16)  in  a  manner  similar  to  that  of  Danbury  Brown  and  Twiss. 

However,  to  develop  an  approach  adaptable  to  arbitrary  orders  of  corre¬ 
lations,  as  wc  as  to  allow  consider;  tion  of  scattering  surfaces  with 
arbitrary  roughness,  we  start  by  examining  a  relationship  valid  for  real 
fourth-order  gaussian  processes’’  (here  in  the  real  electric-field 
spatial  variable),  where 

2  (V'Cxp  V  (x.,))2  =  <AV2(Xj)  AV2(x?)>  (17a) 

=  (AI(Xj)  A I  (xn))  1 17b) 

and  where  we  define 

AV2(xJ  n  V2(x  )  -  <V?(x  )>  (18a) 

— 51  II  “11 

a nc 

AI  ( c  )  =  I  (x  )  -  <  I  (x  )>  .  (ISb) 

-n  -n  -n 
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Equation  (17b)  follows  from  the  definition  of  intensity,  and  the  cngle 
brackets  (>  indicate  an  ensemble  (spatial)  average,  not  the  more  usual 
time  average. 

Using  the  identity 

<AI(xp  AI(x2)>  =<I(xp  I(x2)>  -UUj))  (I(x2)>  ,  (19; 


we  can  write 


<Kx1)  I(x2)> 


=  2(V(x1)  V  (x2)>  +  UQ^))  (I(x2)> 


We  can  conclude  from  Eq.  (20)  that  the  second-order  intensity  correlation 
is  composed  of  two  terms,  of  which  one  forms  the  square  of  the  second- 
order  field  correlation;  the  other  is  a  spatial  dc  term,  of  no  value  here. 
Hanbury  Brown  and  Twiss  eliminated  a  similar  temporal  term  by  means  of  a 
dc  block  in  their  electronic  apparatus.  The  mathematical  origin  of  this 
term  can  be  observed  in  Eq.  (16)  for  the  case  p  =  0,  corresponding  to 
the  infinite  time  average.  For  p  =  0,  the  product  of  intensities  in  the 
far  field  for  just  two  points  is  a  constant;  but,  as  the  mean  position 
of  the  points  is  translated  spatially,  the  product  varies  and  thus  for 
the  spatial  averaging  case,  terms  in  p  =  0  must  be  evaluated. 

Similarly,  for  the  complex  field  variable  we  can  write,  using 
Eqs.  (6)  and  (16), 


UQc^T)  I(x2,T)>  =  |<r(x1,x2;  0;T)>  |‘ 


(21a) 


CO  00 

1  I  2 

—  — r  o  do  dp  sinc(pT)H(p)  exp(ipr/c) 

w  i  L 


||  <A(C1,a)  A(£2,o)  expti^-^p)]) 
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where  we  have  interchanged  the  order  of  integration  and  averaging. 

Next,  we  make  the  following  transformation  to  spatial  center-of- 
mass  coordinates  where 

1  =  ij-I2  £  =  j  .  (22) 

Introducing  these  into  liq.  (21b),  following  some  algebra,  we  find 
(suppressing  constant  terms) 


<I(x1,1)  1(x2,T)>= 


(cr) 


o2do 


sinc(pT)i/(o)dp  exp(ipr/c) 


X 


exp[i<J.(f,p)]> 


*  exp 


[  (A+i/2)  * 


X 


K&-i/2)-x7] 


(23a) 
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In  going  from  Eq.  (2^a)  to  Eq.  23b),  we  have  dropped  the  Fourier- 

transfcrm  terns  in  the  difference  frequency  variable  p,  since  they  are 

clearly  negligible  in  the  far  field.  The  ensemble  average  nas  been 

expressed  as  a  product  of  two  terms:  C(f)  is  a  normalized  correlation 
14 

function  describing  the  coherence  interval  over  a  rough  surface,  and 
I(£,«u)  is  th-  intensity  in  the  global-spatial  variable.  We  have  made 
the  reasonable  assumption  that  the  field  amplitude  is  constant  within 
a  given  coherence  area  of  the  source;  specifically,  for  f_  sufficiently 
small  that  C(_f)  ^  0,  A(g+f/2,c0  =  A(£-f/2,o)  =  A(£,cr).  This  final 

approximation  is  made  under  the  assumption  that  the  coherence  area  of 
the  source  is  small  relative  to  the  total  source  area,  an  assertion 
already  made  in  an  earlier  argument  for  gaussian  statistics.  In 
Eq.  (23c),  the  hat  ~  indicates  a  two-dimensional  spatial-Fourier 
transform. 


In  order  to  describe  a  spatially  incoherent  surface,  C(f)  is 
usually  allowed  to  assume  the  role  of  a  6  function.  Thus,  Eqs,  (23) 
are  reduced  to  a  single  integral  in  two  space.  Given  this  form,  we  see 
from  Eqs.  (23)  that  the  ensemble-averaged,  two-pont  intensity  corre¬ 
lation  in  the  far  field  is  proportional  to  the  modulus  of  the  spatial 
Fourier  transform  across  a  spatially  rough  surface,  assuming  a  suffi¬ 
ciently  short  exposure  time  2T.  However,  only  for  the  case  that  the 
intensity  distribution  over  the  source  has  even  symmetry  :an  the  phase 
of  the  spatial  transform  be  inferred  and  used  to  invert  uniquely  Eqs. (23) 
to  derive  the  intensity  distribution  on  the  source,  I(^,u). 

Finally,  the  form  of  En.  (23c)  shows  explicitly  that  the  two-point 
ensemble-average  intensity  correlation  in  the  far  field  is  proportional 
to  the  product  of  two  spatial  power  spectra.  The  spatial  power  spectrum 


is  multiplied  by 


the  spatial  power  spectrum  of  the  correlation  function  describing  the 
surface  roughness.  If  the  surface  is  sufficiently  rough  that  this 
function  approximates  a  6  function,  then  the  transform  is  essentially 
constant,  and  all  spatial  frequencies  of  the  source  can  be  inferred, 
ricwever,  as  the  correlation  interval  increases, 


acts  to  band  limit  the  detectable  spatial  spectrum  of  the  source.  This 

17 

effect  is  discussed,  for  example,  by  Kinsly  for  the  case  of  micro¬ 
densitometer  imaging  with  partially  coherent  light. 
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V.  AN  EXPERIMENT 


To  illustrate  the  above  ideas,  we  suggest  a  simple  experiment 

embodying  these  mathematical  ideas.  A  helium-neon  laser  in  single- 

axial-mode  configuration  is  used  to  transilluminate  a  symmetrical  source 

(Fig.  1).  The  source  has  a  random-phase  character  to  obtain  spatial 

incoherence.  The  far-zone  intensity  pattern  is  recorded  by  film,  using 

1 8 

an  exposure  time  less  than  the  reciprocal  of  the  source  bandwidth. 

Next,  the  film  is  developed  so  that  it  is  linear  in  intensity  and 
used  to  make  two  identical  positive  transparencies.  The  positives  are 
then  placed  in  a  collimated  beam  (Fig.  2)  to  form  the  correlated  in¬ 
tensity  over  an  averaging  area.  The  signal  transmitted  by  the  trans¬ 
parency  pair  is  optically  Fourier  transformed,  a  dc  block  is  inserted 
to  remove  the  unwanted  average  term,  and  the  total  remaining  irrad’ance 
is  measured.  This  signal  represents  the  mathematical  expression  given 
by  Eqs.  (23)  for  the  transparency  spatial  lag,  d  =  Xj-x,,.  Since  the 
source  is  known,  a  priori,  to  be  symmetrical ,  the  transform  of  the 
source  intensity  is  real.  The  square  root  of  the  correlation  signal  is 
proportional  to  the  spectrum,  which  is  then  known  as  a  function  of 
spatial  la’.  Finally,  this  two-dimensional  signal  is  Fourier  trans¬ 
formed  by  machine  to  give  the  scaled  source  irradiance. 

VI.  SUMMARY  AND  CONCLUSIONS 

Developing  an  intermediate-average,  mutual  coherence  function  as 
a  starting  point,  we  have  derived  an  expression  for  the  two- point,  in¬ 
tensity  ccrrelation  in  the  far  field,  independent  of  time  averaging 
except  for  the  temporal  resolution  of  the  detector.  This  result  is 
valid  for  narrow-band,  high-intensity  light  scattered  from  a  spatially 
rough  surface  of  arbitrary  coherence  area. 

There  are  a  number  of  special  benefits  from  detecting  images  by 

the  technique  of  intensity  correlation,  (i)  The  method  is  relatively 

2 

insensitive  to  the  effects  of  atmospheric  scintillation.  (ii)  Because 
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A  pair  of  transparencies,  T  and  T  ,  are  made  from  film  record  in 
plane  F  of  Fig.  1.  Intensity  autocorrelation  function  is  measured 
by  power  meter  for  spatial  lag  d^-x^;  dc  b1ock  is  inserted  between 
lenses  at  arrow  to  remove  unwanted  signal. 


the  signal  is  detected  in  the  spatial -transform  domain,  high-frequency 
detail  about  the  scattering  surface  translates  to  large  spatial  lags  in 
the  far  field.  This  result  could  be  particularly  important  at  frequen¬ 
cies  where  detector  resolution  is  not  well  developed.  Ciii)  A  special 
advantage  to  intensity  interferometry  in  the  spatial  domain  is  the 
utilization  of  gaussian  statistics  in  the  spatial  (not  temporal)  sense. 

By  this  method,  sources  with  nongaussian  time  statistics  (such  as  single- 
axial-mode  lasers)  can  be  utilized,  (iv)  Still  another  advantage  of 
spatial  detection  is  that  images  of  moving  surfaces  can  be  formed  using 
brief  exposures. 

We  have  therefore  shown  that,  given  a  symmetrical,  spatially  in¬ 
coherent  source  illuminated  by  hijh-intensity  light,  the  far-zone  in¬ 
tensity  pattern  can  be  used  to  form  the  optical  image  of  the  source  if 
the  signal  is  recorded  with  sufficiently  short  time  resolution. 
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